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The dynamics of an ensemble of linear disturbances in boundary-layer flows at various Reynolds numbers 
is studied through an analysis of the transport equations for the mean disturbance kinetic energy and energy 
dissipation rate. Effects of adverse and favorable pressure-gradients on the disturbance dynamics are also 
included in the analysis. Unlike the fully turbulent regime where nonlinear phase scrambling of the fluctuations 
affects the flow field even in proximity to the wall, the early stage transition regime fluctuations studied here 
are influenced across the boundary layer by the solid boundary. The dominating dynamics in the disturbance 
kinetic energy and dissipation rate equations are described. These results are then used to formulate transition- 
sensitized turbulent transport equations, which are solved in a two-step process and applied to zero-pressure- 
gradient flow over a flat plate. Computed results are in good agreement with experimental data. 


I. Introduction 

A current challenge in the prediction of wall-bounded turbulent flow fields within the realm of single-point closures 
is to consistently predict the (upstream) flow field and location of the natural transition region. The most widely used 
method for predicting natural transition is the e A method (see, for example, Drela and Giles 1 ), combined in some 
fashion to a turbulent flow solver. But this technique can be cumbersome and often requires considerable effort on the 
part of the user. 2 

Schmidt and Patankar 3 analyzed the natural capability of low-Reynolds-number K-s models to simulate boundary- 
layer transition. Although there were significant limitations regarding accuracy, the models were considered to be 
useful as predictive engineering tools. Wilcox 4 took advantage of the capability of two-equation models to describe 
nonlinear growth of flow instabilities to develop a modified K-uj model that predicts transition. Walters and Leylek 5 
developed a single-point method that accounted for the growth of Tollmien-Schlichting waves (for natural transition) 
and also employed an additional transport equation to represent the growth of nonturbulent streamwise fluctuations 
(for bypass transition). Wang and Perot 6 applied a non-equilibrium turbulent potential model to both natural and 
bypass transition. 
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Examples of other transition-sensitized turbulence models, applicable to naturally transitioning (lows, exist in the 
literature, 7-10 but their development has been based more on empirical grounds. In contrast, transition-sensitized 
turbulence models, applicable to intermittent Hows 11 and flows with bypass transition, 12 * 13 have been extensively 
studied and refined over the last decade. 

Motivated by the success of models in flows with bypass transition, a methodology is pursued here that can lead 
to the development of models applicable to generic flows in a natural, low-disturbance environment. In order to 
develop a general linkage between the transitioning laminar flow and the turbulent flow in a developing boundary 
layer, it is necessary to have a common mathematical framework from which the disturbances in both regimes can be 
described. In previous studies, 14 * 15 such a framework was developed by coupling a deterministic description of the 
evolution of disturbances in the laminar regime with an analysis of an ensemble of such disturbances. The approach 
was based on the observation that, even in the laminar regime, every flow is subject to an inevitable uncertainty in 
initial conditions. Therefore, although each individual disturbance evolves deterministically, a probability distribution 
must be introduced for the calculation of ensemble mean properties of the various disturbance statistical moments. 

The mathematical methodology developed previously is now extended to a spatially developing boundary layer. 
This study focuses on the earliest stages of transition in a low disturbance environment where the disturbance ampli- 
tudes are small enough that nonlinear interactions can be neglected. Most current Reynolds-averaged Navier-Stokes 
(RANS) type closure models have difficulty representing pre-transition flow fields and transition location, particularly 
with nonzero pressure gradients or for three dimensional mean flows. Although three dimensional mean flows are not 
considered here, the effects of zero, favorable, and adverse pressure gradients on the kinetic energy and dissipation 
rate budgets are considered. The information obtained from these balance equations is used to guide the formulation 
of a unified model, which is then applied to the flow over a flat plate in zero pressure gradient. 


II. Theory 


A flat, solid-walled boundary is considered. The incompressible flow analysis focuses on downstream locations 
away from any leading-edge effects. The three-dimensional, laminar disturbance modes under consideration are 
bounded solutions of the linearized Navier-Stokes equations. Linear disturbance fields are generated that are deviations 
from an ensemble-mean boundary-layer velocity field given by the Falkner-Skan family of solutions, U(Re, 0 ,x 2), 
applicable to flow fields with zero, adverse and favorable pressure gradients. The effects of pressure gradient are in- 
troduced through the parameter 0 : 0 ~ 0, zero pressure gradient; 0 < 0, adverse pressure gradient; 0 > 0, favorable 
pressure gradient. 

The velocity and pressure disturbance fields can be constructed from disturbance mode solutions by considering 
an ensemble of linear superpositions of modes with mode amplitudes <£(u;, fc 3 , xi 0 ) so that 


duj dk^dx 10 fc 3 , a>i 0 )ui(x2,u;, fc 3 , i?e)e 


i( f 1 kidx\ 4- k$xj, — wf) 


«i(x,t) = J 

p(x,t) = J dk 3 dxi 0 $(u},k3,xio)p(x2,u),k3,Re)e { ^w k ' dXl+ kjX> 


(1) 

( 2 ) 


where (k\(Re, w, fc 3 ), 0, fc 3 ) is the wave vector in coordinate directions (xi , X2, x 3 ) associated with the most unstable 
mode, (ki complex and & 3 real), uj is the (real) frequency, xio is the position where the disturbance mode originates, 
and Ui and p are the disturbance velocity and pressure profiles for the most unstable mode. This ensemble is de- 
scribed by a probability distribution for the amplitudes so that their mean, ($(u;, fe 3 , £10)) is zero, and covariance is 
homogeneous (xj and x 3 — directions) and stationary. 


(<I>*(u;, Ar 3 ,£i 0 )$(u/, k' s , x' 10 )) = 6 (u> - u/)£(/c 3 - k^)S(x 10 - x' 10 )P(^, k-$, £10). (3) 

Assuming that fc 3 , and x i0 are independent random variables, the probability density P(u\fc 3 ,x 10) is partitioned 
as 

P{u>,k 3 ,x 10) = Pu,(u;)Pk 3 (k 3 )Px w (xio)- ( 4 ) 
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The frequency space probability density P^( u?) is assumed to correspond to band limited white noise with a flat 
spectrum 


PJlo) 


1 


~ u ; l 

0 


uJL < ^ < &H 
otherwise 


(5) 


where — 0.04 and u )h = 0-08 are chosen to encompass the region of disturbance growth computed from linear 
stability theory. For the spanwise wavenumber probability density, the central limit theorem argues in favor of a 
Gaussian 

jww - 75 ^“ p (-£) <61 

with variance a 2 large enough to encompass a sufficient range of energetic oblique modes. It was found that the 
results did not change for a > 0.1. In this study a ~ 0.1 was chosen and ks values ranged from —0.3 to 0.3. The 
disturbance modes are assumed to originate with equal probability along the boundary layer in accordance with the 
assumed homogeneity of the freestream turbulence field. In order to ensure that the same disturbance ensemble is 
generated at each streamwise location xj, the probability density function P Xl0 (x i 0 ) is chosen to have the form 


Pr lO (*10) = S{X X -X 10 ) 


(7) 


This eliminates the exponential growth factor in Eqs. (1 ) and (2), and ensures that each disturbance mode has the same 
weight. 

The second-moment correlations can be formed from the disturbance velocity and pressure fields given in Eqs. (1) 
and (2). For example, the disturbance stress tensor is given by 

Tij{x i.X 2 ) = u’Uj + UiU* ) 


1 

2 


/ 


(Lj dk 3 P^(m)P k3 (k 3 ) ( u*Uj 


+ Uiiij) , 


with the disturbance kinetic energy, K = tu/ 2, and the disturbance isotropic dissipation rate is given by 


c(x 1 .X 2 ) 


J_ / dttj 
Re \ dx k dx k / 


( 8 ) 


= PjdLxdksP^PkM 


X 



du* duj 1 
3 *** dx 2 dx 2 j 


u,u 


(9) 


Note that, although streamwise homogeneity in the choice of ensemble has been maintained, all disturbance cor- 
relations retain a dependence on x\ through either a direct Reynolds number dependence (as in the disturbance dissi- 
pation rate) or through the implicit Reynolds number dependence of u 7 and k x . Finally, in order to compare results at 
different values of Re$ and 3, the kinetic energy of the disturbance field at each Ree and 3 is normalized such that 
the integral of the kinetic energy across the boundary layer in x 2 is unity. It is now possible, from these quantities, to 
construct a detailed mapping of the ensemble-averaged disturbance field. 


A. Disturbance energy budget 

The derivation of the transport equations for the disturbance kinetic energy, K — r„/2, follows directly from the 
linearized Navier-Stokes equations coupled with an average over an ensemble of disturbances described by the distri- 
bution P( a?, fc 3 , xio). The resulting transport equation can be written as 

DK 

— =P-£-bi7-b£> (10) 
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where 


V = 



( 11 ) 


77 

V 


1 d 

2 dxj 


(UjP + UjP*) 


J_ d?K_ 

Re dxjdxj 


(12) 


(13) 


are, respectively, the disturbance energy production, pressure-transport of disturbance energy, and viscous diffusion of 
disturbance energy. The disturbance energy dissipation rate is e given in Eq. (9). 

Using the solutions u ? and p (Eqs. (1) and (2)) obtained from the linearized disturbance equations coupled with 
the probability density distribution (Eq. (4)), the components of the disturbance energy equation can be obtained. It 
is instructive to examine the behavior of the disturbance ensemble at different streamwise locations. At low Reg (not 
shown) where there are no growing modes in the ensemble, the disturbance kinetic energy decays ( DK/Dt < 0) 
because production, the contribution from dissipation and viscous diffusion are all negative. Only pressure-transport 
contributes positively to the energetic balance over most of the layer. Very near the wall, the dynamic balance is 
between the kinetic energy dissipation rate and viscous diffusion, which is analogous to the fully turbulent case very 
near the wall. 

Figures 1 and 2 show the energy budget across the boundary layer at two different values of Reo. At each value 
of Reg , the distributions of energy production, dissipation rate, pressure-transport, and viscous diffusion are shown. 
Here DK/Dt is calculated as an instantaneous quantity, and not as a finite difference. Figure 1 shows the energy 
balance at Reo = 393; a region of the flow where some modes included in the ensemble are unstable. At this value 
of Reo , there exists regions of the boundary layer where, for example, production is positive and pressure-transport is 
negative although overall DK/Dt < 0. Very near the wall the balance is predominantly between viscous diffusion 
and disturbance dissipation rate. Figure 2 shows the energy balance at Reg = 1040. Increasing the Reynolds number 
causes an increase in the magnitudes of all of the terms in the balance. At larger values of Reg, the qualitative features 
of the energetic balance remain unchanged with only the tails of each distribution decaying more slowly with 
reflecting the downstream boundary-layer growth. At this and larger values of Reg . the kinetic energy is growing with 
DK/Dt > 0 across the boundary layer. Again, the near wall balance is between viscous diffusion and disturbance 
dissipation rate. It is apparent from Fig. 2 that, outside of the very near-wall region, the production and pressure- 
transport have opposite sign and are nearly proportional to one another although viscous diffusion also contributes 
to the balance. It is possible to show that this proportionality between the pressure-transport and energy production 
can be established through an analysis of the Poisson equation for the pressure fluctuations. In addition, the role of 
viscous diffusion in the dynamic balance across the boundary layer differs from the turbulent regime and changes with 
Reynolds number. In the fully turbulent case, the viscous diffusion changes sign with distance from the wall; however, 
for the linear disturbances examined here, the viscous diffusion is positive across the boundary layer at Re = 1040, 
and only changes sign at the lower Reynolds number. 

The effect of imposed pressure gradient on the disturbance kinetic energy budget is shown in Figs. 3 and 4 at Reg as 
1040. In the adverse pressure gradient case (/3 < 0), the dynamic balance between the terms in the transport equation 
remain unchanged; however, there is a broadening of the distribution across the layer and an overall substantial increase 
in the amplitude levels associated with each term. For the favorable pressure gradient case shown in Fig. 4, the dynamic 
balance is significantly altered qualitatively from the zero pressure gradient case. The balance is more consistent with 
the Reg = 393 zero pressure gradient case shown in Fig. 1. At this streamwise station, the ensemble of disturbances 
included both stable and unstable modes, but the overall balance showed a decay of disturbance kinetic energy across 
the boundary layer. 

These results confirm that even though the individual disturbance modes satisfy the linear Navier-Stokes equations, 
an ensemble of such disturbances can display a complex range of energetic dynamics that characterize the early stage 
transition regime of a developing flow. In the next section, the dynamic balance of the disturbance dissipation rate is 
examined. 
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B. Disturbance dissipation rate budget 

The transport equation for the disturbance dissipation rate 5 is also directly derivable from the linearized momentum 
equations. The resulting form is given by 


§? = ?\ + + Pr. -T + n € + v e 

V e 


where 

T>\ = 
P 2 = 

t >3 __ 

r e 

T = 

n e = 

2?, = 


1 dU / du\ du‘> ^ du\ du*, \ 
Re dx 2 \dxj dxj ' dxj dxj / 

duj_du^\ 

Re dx 2 \ 9xi 9x2 9xj 9x2 / 

1 d 2 U / „ 9m i 9wj \ 

/?e <2xj \ 2 dx.2 W ” 9x2 / 

2 / as \ 

i?c 2 ^ j 

1 9 / dp duj dp* \ 

7?e 0Xj \ ^ / 

1 d 2 £ 

Re dxidxi 


(14) 

(15) 

(16) 

(17) 

(18) 

(19) 

( 20 ) 


are, respectively, mixed production, production by mean velocity gradient, gradient production, destruction, pressure- 
transport, and viscous diffusion. (The notation used here is consistent with the notation used by Rodi and Mansour. 16 ) 
Analogous to the treatment of the DK/Dt term, De/Dt is not a finite difference but is computed as an instantaneous 
quantity. 

Figures 5 and 6 show the balance of terms across the boundary layer for the disturbance dissipation rate budget. 
With the exception of the destruction and viscous diffusion of dissipation, which essentially balance out across the 
entire boundaiy layer, all terms are negligible. This behavior suggests that in the absence of a production mechanism 
associated with vortex stretching, a nonlinear effect, direct dissipation of disturbance energy is confined to regions in 
close proximity to the wall. In contrast with the disturbance energy balances, the same (qualitative) dynamic balance 
is achieved for the two values of Res shown and was also found for all values of Res at which these calculations were 
carried out (Re# < 2000) The only change that was found was the broadening of the distribution across the layer with 
increasing Reo. While energetic equilibrium ( DK/Dt ~ 0) only occurred at a single value of Reo, Fig. 5 shows 
that De/Dt is very small and « 0 compared to the dominant terms, destruction and viscous diffusion of dissipation, 
throughout the entire range of Reo values studied. In fact, it was found that De/Dt was always at least three orders 
of magnitude smaller than the dominant terms. 

Although not shown, the dissipation rate dynamics are relatively insensitive to the imposition of pressure gradi- 
ent on the flow. This once again shows that the disturbance dissipation rate is not a (significant) contributor to the 
disturbance dynamics in the early stage transition regime. 


HI. Model 


The Reynolds averaged Navier-Stokes approach for predicting turbulent flow fields continues to be the primary 
engineering tool for calculating entire flow fields at high Reynolds numbers. However, closure models ranging from 
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algebraic to differential stress models usually cannot correctly predict the transition process. Nevertheless, in wall- 
bounded flows for example, the success of a turbulence model is often judged by its ability to predict the developing 
flow field. Such an assessment is often clouded by the fact that the region of transition is not well-determined by the 
turbulence model. 

In the previous section, a set of transport equations for the disturbance kinetic energy and (isotropic) dissipation 
rate were analyzed. The disturbance variables were ensemble averages of deterministic disturbances obtained from 
the linear Navier-Stokes equations, and associated with a chosen probability density function, Eq. (4). While the 
probability distribution associated with the linear disturbances may be different than that for the turbulent fluctuations, 
the means associated with both ensembles are formed similarly. Thus, as was done in the last section, quantities such 
as the kinetic energy and dissipation rate and their respective transport equations can be analyzed analogous to those 
associated with the turbulence fluctuations. 

As assumed throughout this analysis, the ensemble of disturbances examined are representative of the disturbance 
field in the early stages of a naturally transitioning wall-bounded flow. 17 The budgets discussed earlier can then be 
used, as in the development of RANS models for turbulent flow, to develop a RANS-type model for the early stage 
transition regime. The low Reynolds number (Reo = 393) results, where the kinetic energy production is negative, are 
not of interest here, since they represent a regime in which disturbances decay instead of transitioning into turbulence. 
Rather, the results for Re& = 1040 with zero and adverse pressure gradient are relevant and will be used. 

The focus of this analysis is on linear eddy viscosity two-equation models. The eddy viscosity is ordinarily 
associated with turbulence, but in the current analysis it is also considered a disturbance viscosity that applies to linear 
disturbances in the early stages of transition. As in the turbulent case, the task then is to determine closure models for 
the destruction, pressure transport, and diffusion terms that appear in the disturbance kinetic energy and dissipation 
rate equations. At the outset, it is actually necessary to only consider models for the pressure- transport and diffusion 
terms. The destruction of disturbance kinetic energy is governed by the disturbance dissipation rate equation itself. 

The base K-e turbulence equations chosen for the purposes of the current investigation are written: 


DK ^ 1 d 

Dt ^ P 9xj 


M -I 

<*K 


OK 

dxj 


= v-c —t + 12L 

Dt K° ElV e2 K 2 + p dxj 


»t\ de_ 

o e )dx n 


where V = — t, :i dU t /c)xj « 'hy-j-SjjSyj is the production term, and the resulting eddy viscosity is given by: 


( 21 ) 

( 22 ) 


HT = pv T = C^fnp— (23) 

The constants are: C e i = 1.44, C s 2 = 1.83, Ok = 1.0, cr £ = k 2 j[^JC^{C £ 2 — C s 1 )], = 0.09, and « = 0.41. The 

function is given by: 


/ 2 = l-e X p(-i§) <24) 

where Rex = py/Kd/p, and d is the distance to the nearest wall. For simplicity, the function (which is typically 
used to insure proper asymptotic behavior of the turbulent variables as the wall is approached) is currently taken to be 
unity. 

Inspection of Fig. 2 shows that the plots of disturbance energy production, diffusion and pressure-transport assume 
similar shapes from £3 ~ 10 to #2 ~ 00, and combine to produce a positive DK/Dt. This suggests that the sum 
of these three terms may be modeled as an effective disturbance energy production Vd • Regarding the dissipation rate 
equation, on the other hand, the data indicate that De/Dt « 0 throughout the layer. These ideas lead to the following 
possible forms for the modeled transport equations for the disturbance energy and dissipation rate: 


DK 

Dt 


r d -e 


(25) 
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£i 

Dt 


0, 


(26) 




where V d represents the sum V -f i7 -h V appearing in Eq. (10) and is modeled as CaV* . For this two equation 
formulation, a Boussinesq type relation is assumed for the disturbance shear stress so that the production V f can be 
written in the form 2 is d SjjSij, with v d the disturbance mode viscosity (y d = C tld K 2 /c). 

In order to bridge together the two flow regimes (represented by Eqs. (21 ) & (22) and Eqs. (25) & (26)), an indicator 
function X is used. X represents the probability that the flow at any given point is turbulent, with 0 < X < 1. Thus, 
X — 0 represents purely laminar flow and X = 1 represents purely turbulent flow. The precise form that the X-equation 
needs to take is an ongoing subject of active research. Currently, a provisional measure has been taken in an effort 
to reproduce the known dependence of transition on freestream turbulence intensity. This form of the X-equation is 
derived below. 

The original K and £ equations are modified using the function X such that when X — 1 the fully turbulent 
equations are recovered, and when X = 0 the laminar disturbance equations are recovered. (However, as an aid to 
improved numerical behavior of the solution, the disturbance form of the equations are modified to retain dissipation 
terms. The destruction term in the dissipation rate equation is also retained in order to prevent too-rapid freestream 
eddy viscosity decay.) The final forms of the K and £ equations are: 


DK 

Dt 

De 

~Dt 


= (1 -X)V d + XV-s + 


U+i— 

LV ffAy 

<)K 
& x j . 

[H?) 

) — 


(27) 

(28) 


Here, V d « 2v* r S ij S ij (C d C» d /C fl ), C d = 0.3, and C» d = 0.026. 

The indicator function X should obey an evolution equation that produces the correct dependence of transition 
location on the freestream turbulence level Tu. As freestream turbulence increases the transition location for a flat plate 
boundary layer is known 18 to move forward along the boundary layer, i.e., to smaller values of streamwise coordinate. 
The indicator function should also be related to some local quantity that gives a measure of the level of turbulence 
within the boundary layer. Motivated by these considerations, the behavior of the turbulence Reynolds number Re t — 
pK 2 / pe was studied for laminar (X = 0) solutions of Eqs. (27) and (28). (Note that Ret is proportional to the 
eddy viscosity pr appearing in the production and transport terms of the disturbance kinetic energy and dissipation 
rate equations.) Figure 7 shows contour plots of Re t for three different levels of T u. The higher the freestream 
Tu , the higher the maximum Re t levels, and the further forward the location where Re t first exceeds the freestream 
level. This suggests that the rate of increase in the indicator function could be related to Re t . As a first provisional 
attempt at constructing the X-equation, a source term proportional to the local value of Re t minus its freestream value 
was introduced. With this choice of source term, X begins to increase when Re t exceeds its freestream value. A 
diffusion term was also included. This equation produced the desired dependency of transition location to freestream 
Tu, but it should be considered work-in-progress. (For example, a destruction term may also be desirable to provide 
relaminarization capability.) The current form of the indicator function transport equation is given by: 


Dt p ox j 


Pt dX 
ox dxj 


(29) 


where V x = Ret — Ret, oc , Cj = 4 and o x = 1. The production term is limited via V x = mmfmax^z, 0), 50], and 
X is limited to be between 0 and 1. Eqs. 27. 28, and 29 are solved in conjunction with the RANS equations. A specific 
procedure is necessary for obtaining consistent solutions to the equations: 


1. First, solve only RANS plus the two disturbance equations (Eqs. 27 and 28 with X set to zero). This step 
establishes the “disturbance levels” of K and £ in the vicinity of the body. It is important that this step be fully 
converged (the final answer is strongly dependent on the level of convergence of this step). 

2. Restart from the solution of the previous step using RANS plus all three equations, Eqs. 27, and 28, and 29. 
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The boundary conditions at solid walls are: K w — 0, (de/dn) w = 0, and l w — 0. At freestream inflow the three 
quantities are taken to be their freestream values, and at freestream outflow they are extrapolated from the interior of 
the domain. 

The freestream conditions are established as follows. For the indi cator function, = 0. Given a freestream 
turbulence intensity Tu (in percent, defined as Tu = 100^2AY(3L7b))> the freestream nondimensional A7o is given 
by: = Koo/a ^ = 1.5[(Tu)(i\/ oo )/100] 2 . In the current methodology, the freestream i?e t>oc is taken to be 

proportional to Therefore, freestream £qo is proportional to K 0 c . One advantage to using this relationship 

is that the resulting decay rate of the turbulence quantities in the freestream is identical regardless of freestream Tu. 
Currently: e ^ = ^oo^oo/{Poo a to) — (5.4x 10 _8 )A^. Table 1 gives freestream levels for various Tu , for = 0.2. 

Table 1. Example freestream levels for M ^ =0.2 


Tu , percent 

R og Too f Q*o o 

= ^ooMoo/(Pooa^o) 

R^t, oc 

/A, OO 

0.042 

1.058 x 10“ 8 * 

5.715 x 10- 16 

0.20 

0.018 

0.1 

6.000 x 10" 8 

3.240 x 10- 15 

1.11 

0.10 

0.2 

2.400 x 10" 7 

1.296 x 10~ 14 

4.44 

0.40 

0.26 

4.056 x 10- 7 

2.190 x 10- 14 

7.51 

0.68 

0.34 

6.936 x 10“ 7 

3.745 x 10- 14 

12.8 

1.15 

0.70 

2.940 x 10" 6 

1.588 x 10- 13 

54.4 

4.90 

1.25 

9.375 x 10“ 6 

5.063 x lO" 13 

173.6 

15.6 


Like the original k-s model itself, the current method does not maintain freestream turbulence levels (i.e., there 
is decay of the turbulence quantities with distance from the inflow boundary, where freestream levels have been set). 
The current model also shares the characteristic that for high freestream Tu (greater than 0.5% or so), the freestream 
Re t , og (or levels can get very large. See Table 1. With unrealistically high background eddy viscosity levels, 
laminar regions in the flow are corrupted. Unfortunately, if one tries to correct this behavior by increasing then 
freestream decay rate becomes much too rapid. The current levels were chosen as a compromise between rapidity of 
decay (current freestream eddy viscosity decays to 74% of its inflow value by the time it reaches the plate leading 
edge) and having excessively high freestream eddy viscosity levels. In any case, the primary focus of this method is on 
natural transition, which should be applicable only for low Freestream Tu. A focus of future work will be developing 
the capability to maintain freestream turbulence levels in the method with little or no decay. 

IV. Results 

The method described in the last section was implemented into the finite-volume RANS code CFL3D of Krist 
et al. 19 Results were generated for a flat plate in zero pressure gradient, at Re = 6 million based on plate length 
and M — 0.2. All results were two-dimensional, described from this point forward by the ( x , y ) coordinate system 
notation. A fine grid of size (385 x 129) was used, along with a medium (193 x 65) and coarse (97 x 33) level created 
by using every other grid point from the next finer grid. The grid extended from an x-location of —1/3 to 1, with the 
plate itself located between 0 and 1. The grid extent in the ^/-direction was approximately 1. The minimum spacing at 
the wall yielded an average minimum y+ of approximately 0.2 on the fine grid, 0.5 on the medium grid, and 1.1 on 
the coarse grid. 

Figure 8 shows skin friction coefficient results plotted as a function of Re x (= u^x/u^q) using the medium grid 
for a range of different Tu , from 0.042% to 1.25%. Transition start locations were taken from the experimental results 
of Schubauer and Skramstad 20 (for 0.042 < Tu < 0.344), and from Abu-Ghannam and Shaw 18 (for Tu — 1.25%). 
Computed results agreed well with the experiments at the intermediate levels (0.1 < Tu < 0.34). Although the 
results at Tu = 1.25% predicted the transition start location reasonably well, the computed Cf levels in general were 
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poor because of the corrupting influence of the too-high freestream eddy viscosity, as discussed in the last section. At 
Tu — 0.042%. the computations predicted transition to be further aft than in the experiment. However, experiments 
often show widely different transition locations at very low Tv levels, due to different noise levels in the wind tunnels. 
See Wells. 21 

Sensitivity of the results to grid size are shown in Fig. 9, for the case with Tu = 0.2%. The transition location 
changed by about 3.5% of the plate length between the medium and fine grids, and by about 5% between the coarse 
and medium grids. These levels of grid dependence are not surprising if one considers that the disturbance equations 
are likely to be quite sensitive to the streamwise gradients. The Re x at the starting location of transition are plotted as 
a function of Tu in Fig. 10; the figure includes the effect of grid density. The experimental “best-fit" transition start 
location curve from Abu-Ghannam and Shaw 18 is also plotted. There was an excellent correspondence between the 
CFD results and the experimental correlation. 


V. Conclusions 

The dynamics of an ensemble of disturbances in boundary- layer flow with and without pressure gradient was 
studied. The individual disturbance modes w ere solutions of the linearized Navier-Stokes equations and the probability 
density function associated with the ensemble accounted for distributions in spanwise wavenumber, frequency and 
initial spatial position. Such an ensemble represents the early stages of a transition disturbance field which ultimately 
leads to a fully turbulent field. 

An analysis of the dynamic balance in the disturbance kinetic energy and dissipation rate equations showed be- 
havior that contrasted with the fully turbulent case at both low Reynolds numbers and in near-wall regions. This new 
insight was critical in properly representing the dynamics in the linear disturbance region. While no self-similar behav- 
ior was reached, the qualitative features of the dynamic balance became independent of Reg and showed that energy 
production, pressure-transport, and viscous diffusion dominated across the boundary layer with a small contribution 
from dissipation rate. The dissipation rate only played a role in very close proximity to the solid boundary where it 
balanced the energy diffusion. 

The results obtained here provided new insight into the dynamic balance of an ensemble of disturbances in the 
early stages of a (naturally) transitioning boundary-layer flow in a low disturbance environment. The information was 
used to develop transport equations for the ensemble of disturbances, which were subsequently coupled using an inter- 
mittency function to corresponding transport equations valid in the fully turbulent regime. The resulting three-equation 
model was solved via a two-step procedure. Results for a zero-pressure-gradient flat plate agreed well with experimen- 
tal data for freestream turbulence intensity levels at and below 0.34%. Transition locations were also predicted in good 
agreement with experiment at higher freestream intensities, but issues related to freestream turbulence decay corrupted 
the quality of those results. Developed from first principles, this method holds promise as a transition-prediction tool 
for the Reynolds-averaged Navier-Stokes equations, with no need for empirical correlations. 
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Figure 1. Disturbance energy balance at different downstream positions (zero pressure gradient), Re$ — 393. 
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Figure 2. Disturbance energy balance at different downstream positions (zero pressure gradient), Reg = 1040. 
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Figure 3. Effect of adverse pressure gradient (J3 = —0.15) on disturbance kinetic energy balance, Re$ = 1045. 



Figure 4. Effect of favorable pressure gradient (/? = 0.15) on disturbance kinetic energy balance, Re$ — 1041. 
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Figure 5. Disturbance dissipation rate balance at different downstream positions (zero pressure gradient), Re# 
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Figure 6. Disturbance dissipation rate balance at different downstream positions (zero pressure gradient), Re e 
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Figure 7. Contour plots of Ret over flat plate (M = 0.2, Re = 6 x 10 6 ) for three different freestream Tu levels, using disturbance 
equations. 
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Figure 10. Re x at transition start location for flat plate as a function of Tu. 
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